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H.E.S.S. observations of massive stellar clusters 
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Abstract. Stellar clusters are potential acceleration sites of very- high-energy 
(VHE, E> lOOGeV) particles since they host supernova remnants (SNRs) and 
pulsar wind nebulae (PWNe). Additionally, in stellar clusters, particles can 
also be accelerated e.g. at the boundaries of wind-blown bubbles, in colliding 
wind zones in massive binary systems or in the framework of collective wind 
or wind/supernova(SN) ejecta scenarios. Motivated by the detection of VHE 
7-ray emission towards Westerlund 2 and assuming similar particle acceleration 
mechanisms at work, Westerlund 1 is an even more promising target for VHE 7- 
ray observations given that massive star content and distance are more favorable 
for detectable VHE 7-ray emission compared to Westerlund 2. Here, H.E.S.S. 
observations of massive stellar clusters in general with special emphasis on the 
most massive stellar cluster in the galaxy, Westerlund 1 are summarized. 



Introduction 



In recent years the field of ground-based VHE 7-ray astronomy experienced a 
scientific breakthrough due to developments in instrumentation, operation and 
data analysis of the 3 rd generation of Imaging Atmospheric Cheren kov Tele- 
scopes (IACTs). T elescope syst ems like H.E.S.S. (Jfflnton et al.l l2004lh MAGIC 



(lLorenz et al.ll2004l V VERITAS (jWeekes et al.ll2002ft or CANGAROO-IH (|Kubo et al 



20041 1 opened a previously inaccessible window for the study of astrophysical ob- 



jects at very high energies. Especially the increase of the number of detected 
VHE 7-ray sources by one order of magnitude is a merit of the Galactic Plane 
Scan (GPS) performed by the H.E.S.S. Collaboration between 2004 and 2008 
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(lAharonian et al.l l2006al : iHoppe et aDl2007l : IChaves et al.l l2008h . In the GPS, 



which covers basically the whole inner galaxy between Galactic longitude 1 
275° and ~ 60°, a rich diversity of astrophysical objects which emit VHE 7-rays 
was discovered. Even if a large part of the source population lacks of a firm 
detection in other wavelength bands, a significant fraction of identified objects 
are connected to the late phases of stellar evolution like e.g. supernova remnants 
(SNRs) or pulsar wind nebulae (PWN). Moreover, one class of objects emitting 
VHE 7-rays could be associated to regions of massive star formation and to the 
birthplaces of the massive progenitors of some SNRs and PWNe, massive stellar 
clusters. 

In this work H.E.S.S. observations of massive stellar clusters with special 
emphasize on the most massive stellar cluster in our galaxy, Westerlund 1 are 
presented. After an introduction of the acceleration mechanisms which are at 
work in such systems (Chapter [2]]), recent results of multiwavelength observa- 
tions and the possible connection of the unidentified source HESS J1614-518 to 
the stellar cluster Pismis 22 are presented in Chapter [3] The detection of Wester- 
lund 2 (HESS J1023-575) in VHE 7-rays is discussed in Chapter H] Finally, the 
detection of VHE 7-ray emission from the vicinity of the massive stellar cluster 
Westerlund 1 is discussed in Chapter [5H 



Particle acceleration mechanisms in massive stellar clusters 



It is widely believed that massive stars form in groups fro m the collapse of ga s 
condensations in clumps, i nside giant molecular clouds (e.g. Hever et al.l (fl 998 ) ; 
Williams. Blitz fc McKeel (feOOOT ) : IPame et al.l (|200ll ): I Jackson et a.1.1 (|200d ^. The 
accretion of cloud material onto the massive stars is stopped, when their winds, 
outflows and UV radiation leads to the dissipation or disruption of the natal 
molecular cloud. Depending on the total mass of the system, these groups of 
stars end as loosely bound associations or as dense gravitationally bound stellar 
clusters. 



2.1. Massive binary systems 

Since a large fraction of massive stars occur in binary (or even triple, quadru- 
ple, etc.) systems (e.g. ( Zinneckerl 120031 : iGiesI 12008) and references therein), 
the collision of their strong superso nic winds with terminal velocities v„o > 



1000 — 5000 km s 1 ( CassineilH 1979), produces strong shocks. Through first- 



order Fermi acceleration el ec trons and protons can be accelerated to high ener - 
gies (lEichler fc Usovi Il993t > ■ iMiicke fc Pohl (|2002l ). iBenaglia fc Romero! (|2003l 1 
and lReimer et al.l (|2006l ) showed, that detectable 7-radiation up to GeV energies 
can be generated through inverse Compton scattering of relativistic electrons in 
the dense photospheric stellar radiation field in the wind-wind collision zone. 
The dete ction of non-ther mal X-ray emission in massive binary systems like e.g. 
WR140 (iDoughertvi 120051 1 has proven that indeed electrons are accelerated in 
the wind collision region to rel ativistic energies . However, just in one hadronic 
emission scenario (proposed by iBednarek ( 20051 )). a detectable signal by current 
and future IACTs is predicted. Here, the inelastic scattering of relativistic nu- 
cleons on particles in the dense stellar wind produces a significant amount of 
neutral pions which then decay into VHE 7-rays. Though, at these high ener- 
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gies, 7-rays suffer from 77 absorption which will diminish the observable flux 
from a close binary system. 



2.2. Collective stellar winds 



Given the fact, that massive stars mostly occur in associations and stellar clus- 
ters, the interaction of their strong supersonic winds leads to the creation of a 
large collective bubble, also called sup er bubble (SB), which is filled with a hot 



and tenuous plasma (e.g. ( Weaver et al.l fl977l ) ) . Turbulent particle accelera- 



tion can then take place, where the interactio n of wind material forms regions 
of strong turbulence and MHD waves (e.g. (jBvkov 120011 : iParizot et al.1 12003 : 

it 



Higdon k Lingenfelter|[2005l : lDwarkadasj l2008^. 



2.3. SN explosions in stellar clusters 

Since massive stars (M > 8Mq) end their lives as SN explosions after a few Myrs, 
an additional contribution of kinetic energy is available to accelerate particles to 
very high energies. The SNR shell will grow quicker due to the lower density in 
a SB and in a medium with higher sound speed, due to the higher temperature. 
This may result in effici ent particle acceleration at MHD turbulences and at 
the b oundary of the SB (jParizot et al.ll2004l : iTang & Wangll2005l : iFerrand et af 
20091 ). 



2.4. Implications 

In a simplified picture, the aforementioned acceleration mechanisms infer basic 
properties of stellar clusters, which in principle could indicate the fraction of 
total kinetic energy in non-thermal particles and therefore a possible emission 
in VHE 7-rays. One obvious property is the total mass of the cluster, since 
this directly determines the number of massive stars and therefore the total 
available mechanical energy in stellar winds (and later in the expanding shells 
of SN explosions). The age of the stellar cluster (and its mass) further implies 
whether or not the most massive members already evolved into SNe and finally, 
the fraction of massive stars located in binary systems defines the contribution 
from colliding wind binaries (CWBs). 

In the next Chapters H.E.S.S. detections of VHE 7-ray sources which could 
be associated to stellar clusters of different mass, age and distance will be pre- 
sented. 



3. HESS J1614 518 



HESS J1614-518 is one of the brightest (~ 25% Crab flux, E > 200 GeV) uniden- 



tificd VHE 7-ray sources discovered during the Galactic plane scan (jAharonian et al. 



2006a). The region was observed for a total of ~ 13 hours resulting in the detec- 



tion of an extended, elliptical shaped emission with a half-width of ~ 0.4° along 
the semimajor axis (Fig. (TJ left). The overall morphology of the source is well 
described by a double peak Gaussian profile (x 2 = 0.7). The differential energy 
spectrum is compatible with a power law (dN/dE= $0 • (E/lTeV)~ r ) with a 
photon index of T = 2.26 ± 0.05 s t a t ± 0.06 sys t and a normalization at 1 TeV of 

\yst) 



$0 = (7.83 ± 0.40 stat ± 0.80 syst ) x 10~ 12 TeV" 1 cm" 2 s" 
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An analysis of the MSX IR data shows several HII re gions and molecular 



cloud complexes flanking the observed VHE 7-ray emissio n (jAmaral et al.lll991 



Matsunaga et al.ll200ll ; iRusseil et al.ll2005l ; iLi et al.ll2006l ). A search for poten- 



tial counterparts reveals several interesting candidates (see Fig. [TJ right). On 
the basis of spin-down power of the pulsars in the field-of-view of HESS J1614- 
518 none of these is powerful enough to account for the total observed VHE 
7-ray emission. Nevertheless, a minor contribution is not excluded. Of the 
Wolf-Rayet (WR) star WR74, the Ant Nebula and the open stellar cluster 
Pismis 22, the latter, located at a distance of d~l k pc, is the most pro mising 



candidate to explain the observed H.E.S.S. emission (|Rowell et al.ll200Sft . Un 



der the assumpti on that 10 stars of B -type (with 20% efficiency) are associated 
to Pismis 22 (see Rowell et al . (2008)), their stellar wind luminosities would be 
sufficient to power the o bserved VHE 7-ra y emission. Additionally, with an es- 
timated age of 40 Myr ( Piatti et alj|2000l ). the system is old enough, that the 
most massive member stars could have already evolved into SNe, giving an ex- 
tra contributio n to the energy ava ilable for part icle acceleration. Furth ermore, 
XMM-Newton (|Rowell et al.ll2008h and Suzaku (|Matsumoto et al.1 120081 ) obser- 
vations revealed an extended non-thermal X-ray emission towards the northern 
part of the source. Remarkably, the Suzaku 'Src B' also shows an extended 
non-thermal nature and is centred on Pismis 22, si milar to other open c lus- 
ters connected to VHE 7-ray emission , like Cyg-OB2 ( Aharonian et all 20051 ) or 
Westerlund 2 (| Aharonian et al.1 [20071 ). 



4. Westerlund 2/HESS J1023 575 



The prominent giant H II region RCW49 and its ionizing cluster Westerlund 2 
are located towards the outer edge of the Carin a arm in our Milky Way at 
a distance of 8 kpc, following iRauw et al.1 ( 20071 ) . Recently Chandra d iscov- 
ered ~500 point sources in the vicinity of RCW49 (|Tsuiimoto et al.ll2004f ). with 
~100 o f them spatially coinc ident with the central open stellar cluster Wester- 
lund 2 ( Townslev et al. 2004! ) . Mid- infrared meas urements with Spitze r revealed 
still ongoing massive star formation in RCW49 (jWhitnev et al.l 12004). The re- 
gions surrounding Westerlund 2 appear evacuated by stellar winds and radiation, 
in contrast to the dust distribution in fine filaments, kn ots, pillars, bubbles, and 
bow shocks throughout the rest of the H II complex (jChurchwell et al.l 12004 ; 
Conti fe Crowther 20041 ). Radio continuum observations by ATCA at 1.38 and 
2.38 GHz indicate two wind-blown shells in the core of RCW49 (jWhiteoak fc Uchidal 
1997), one surrounding the central stellar cluster, the other a Wolf-Rayet star 
WR20b. Westerlund 2, is centrally located within RCW49. It contains an ex- 
traordinary ensemble of hot and massive stars, presumably at least a dozen OB 
stars, and two remarkable Wolf-R ayet stars. The bin ary c haracter of WR20a wa s 
only recently established. Both IRauw et al.l (|2004l ) and iBonanos et all ( 20041 ) 
presented solutions for a circular orbit with a period of ~3.7 days. The derived 
inclination angle of (74.5 ± 2.0)° implies masses of 83.0 ±5.0 and 82.0 ± 5.OM 
for the primary and secondary component, respectively ( Rauw et al.ll2004l ). This 



These stars are in late phases of the evolution, exhibiting strong winds and large mass-loss 
rates. 
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Figure 1. Left: Exposure corrected skymap of the VHE 7-ray excess, 
smoothed with a Gaussian (a = 4.2'). The green contours represent signifi- 
cance levels of 3, 5, 7, 8, and 9cx, after integrating events within an oversam- 
pling radius of 8 =0.1°. The bottom left inlay shows how a point-like source 
would have been seen by H.E.S.S. for this analysis after Gaussian smoothing. 
Right: Finder chart for various objects discussed in the text. Black unfilled 
squares depict X-ray sources. Member and surrounding f ield s tars are also 
shown for the open cluster P ismis 22 (from Piatti et al.l (|2000l ). labeled as 
Pismis 22). Figure taken from lRowell et al.l (|2008h . 



puts the Wolf-Rayet binary WR20a as the most massive of all measured binary 
systems in our Galaxy. 

H.E.S.S. observed the region around Westerlund 2 in 2006 where a point 
source analysis on the nominal position of WR 20a resulted in a clear signal with 
a significance of 6.80". Further investigation s revealed an extended excess with a 
peak significance exceeding 9a (Fig. [2]reft, ( Aharonian et al. 2007)). The source 
is clearly extended beyond the nominal extension of the point spread function 
(PSF), a fit of a Gaussian fold ed with the PSF of the H.E.S.S. instrumen t gives 
an extension of 0.18° ± 0.02° (| Aharonian et al.ll2007l : iReimer et a"fll2007l ). The 
differential energy spectrum for VHE 7-rays inside the 85% containment radius 
of 0.39° can be described by a power law (dN/dE= <J> ■ (E/lTeV)~ r ) with a 
photon index of T = 2.53 ± 0.16 s t a t i 0.1 sys t and a normalization at 1 TeV of 
$0 = (4.50±0.56 sta t ±0.90 syst ) x 10~ 12 TeV" 1 cm" 2 s~ x . The integral flux above 
380 GeV is (1.3 ± 0.3) xl0- n cm- 2 s _1 . 

The detection of VHE 7-ray emission from the vicinity of Westerlund 2 
can be explained by different acceleration mechanisms. One possibility is the 
conversion of kinetic energy available in the colliding-wind region of the binary 
system WR20a into the acceleration of particles to very high energies. Given the 
fact, that no significant flux variability could be detected in the data set and since 
the observed excess is extended with respect to the H.E.S.S. PSF, an association 
between WR20a and the VHE 7-ray emission is not striking. Alternatively, the 
emission could arise from the collective stellar wind effects in the Westerlund 
2 cluster. As discussed before, particles are accelerated by multiple shocks and 
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Figure 2. Left: H.E.S.S. 7-ray count map (oversampling radius of 0.12°) 
of the Westerlund 2 region. The inlay in the lower left corner shows how 
a point-like source would have been seen by H.E.S.S. WR 20a and WR 20b 
are marked as filled triangles, and the stellar cluster Westerlund 2 is repre- 
sented by a dashed circle. Right: Significance contours of the 7-ray source 
HESS J1023-575 (5, 7 and 9cr ), overlaid on a MOST radio image. The wind- 
blown bubble around WR 20a, and the blister to the west of it can be seen as 
de pressions in the radio conti nuum. The blister is indicated by white dots as 
in lWhiteoak fc Uchidal (|1997f ). and appears to be co mpatible in di r ection and 
location with HESS J1023-575. (Figure taken from lReimer et all (|2007f )V 



interact with the dense photon field of the stellar winds. However, the VHE 
7-ray spectrum is expected to be similar that of SNRs at energies below 1 TeV. 
Particle acceleration by MHD motions (shocks, turbulence, etc.) of magnetized 
plasma produced by supersonic flows, which th en penetrate into a dense medium 
may also be important under these conditions (Bykov 200]]). The massive stellar 
winds of Westerlund 2 could ensure sufficient particle injection into the turbulent 
plasma, feeding of magnetic turbulence with energy by wind-wind interactions 
of the massive star association, and the allocation of enough target material for 
the VHE 7-ray production. 



5. Westerlund 1 

Motivated by the detection of VHE 7-ray emission towards Westerlund 2 and 
assuming similar particle acceleration mechanisms at work, Westerlund 1 (Wd 1) 
is an even more promising target for VHE 7-ray observations given that massive 
star content and distance are more favorable for detectable VHE 7-ray emission 
compared to Westerlund 2. 

Westerlund 1 is the only known super star cluster in our galaxy. With 
a total mass of ~ 6 x 10 4 Mq, it is currently the record holder in terms of 
its rich content of stars in the Wolf-Rayet phase. At least 24 WR stars ar e 
known of which >70% are expected to be in binary systems ( Groh et al.| [2006) . 



Additionally the existence of more than 80 blue super-giants, 3 red super-giants, 
one luminous blue variable and 6 (out of 12 in the whole galaxy) yellow hyper- 
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giants is reported in several papers (summarized in Muno et al. ( 2006bl )). The 



dissipated power in the form of kinetic energy in the wind of the WR star s 
alone is L\y ~ 10 39 ergs s" 1 . Given an age of ~ 5 Myrs (|Crowther et al.1 l200d ) . 
the most massive stars in Wd 1 evolved into SNe. Assuming furthermore a SN 
rate of 10 -4 yr -1 , the available kinetic energy for particle acceleration would be 
quadrupled ( Muno et al.ll2006bl ). 

X-ray observations with the Chandra satellite have revealed a magnetar can- 
didate (AXP CXOU J16 47-4552, P = 10.6 s) with a massive progenitor which 
is associated with Wd 1 ( Muno et al.ll2006al ). A detailed analysis also revealed 
extended emission (of the order of arc minutes) which deviates from the typical 
thermal emission that can be seen from many other stellar associations. With 
increasing distance to the cluster core, the hard non-thermal emission dominates 
and line-emission disappear s. The total lumino sity observed in X-rays amounts 
to L x = 3 x 10 34 ergs s~ 1 ( Muno et al.l l2006bl ). representing just 10~ 5 of the 
total mechanical luminosity in this system. The same authors discuss various 
possibilities to explain the 'missing energy' in this powerful system, amongst 
others, the dissipation beyond the Chandra field-of-view by a large scale outflow 
and radiation at other wavelengths. Recent measurem ents suggest a d i stance 
of t he cluster of 4- 5 kpc . Whereas the estimates by iBrandner et al.l ( 20051 ) 
and ICrowther et all ( 20061 ) depend on the star c o ntent of the cluster and there- 
fore suffer from extinction, iKothes Sz Dougherty! ( 20071 ) investigated the neutral 
environment of Wd 1 to find emission and absorption features in the HI data 
which could be connected to the stellar cluster. Interestingly they find two small 
expanding bubbles with dynamical ages of ~0.6 Myrs and a large interstellar 
bubble with a size of 50 pc and a dynamical age of 2.5 Myrs at a distance of 
3.9 ±0.7 kpc. 

H.E.S.S. observed the Westerlund 1 region from 2004 - 2007 during the GPS 
and in pointed observations for 14 hrs under good weather conditions (according 
to the Standard quality selection ( Aharonian et af1l2006bl )) with the full array 
at zenith angles below 55°. Another 22 hrs of good quality data was obtained in 
observations conducted between May - August 2008, resulting in a total data set 
of 34 hrs live-time 0. Preliminary and previously unpublished H.E.S.S. results 
of these observations reveal a spatially extended emission region of VHE 7-rays 
with a total significance of > 15cr, within an integration radius of 1° around the 
Wd 1 position (85% containment) (Fig. [3j left). 

The extension of the VHE 7-ray emission exceeds 2° in diameter, making it 
one of the largest structures observed in VHE 7-rays so far. A total of > 2300 
7-rays are detected within the 34 hours of livetime. The threshold of the analysis 
is at 680 GeV. 

Among the 'established' classes of counterparts for VHE sources like SNRs 
and PWNe only PSR J1648-4611 is coincident with the observed 7-ray emission. 
Given its E/d 2 of 6.2 x 10 33 ergs s _1 kpc -2 it can contribute to the VHE 7- 
ray emission. Recently, the LAT instrument on board of the Fermi satellite 
detected unp ulsed emission from a region coincident with the pulsar position 
(|Abdo et ali l2009'). The Low Mass X-ray Binary (LMXB) GX 340+0 and other 
point-like counterparts are unlikely, given the rather extended nature of the 



2 observation time, corrected for the dead-time of the system 
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Figure 3. Left: Preliminary gaussian smoothed H.E.S.S. 7-ray count map 
(er = 0.13°) of the region around Wd 1 (indicated as white star), filled green 
and black triangles mark INTEGRAL sources and the pulsar PSR J1648- 
4611, re spectively. C yan circles represent SNR candidates from Greens cat- 
alogue (|Greenl [20041 ) . The white circ le indicates the Fermi bright source 
OFGL 1648.1-4606 (|Abdo et al.ll2009ft . The bright region on the right is 
the known H.E.S.S. source HESS J1640-465. Right: HI channel map of 
the Southern Galactic Plane Survey (SGPS) at a velocity of -55 km s _1 
(jMcClure-Grifnths et al.l 20051) corres ponding to a distance of ^4 kpc as ob- 
tained bv lKothes fc Dougherty! (|2007h . Color scale from 20 K to 100 K, over- 
laid are the preliminary H.E.S.S. significance contours from 4 to 8 in steps 
of ler after integrating events within 0.22° radius (as done for the analysis of 
slightly extended sources compared to the H.E.S.S. PSF in the GPS). 



VHE 7-ray emission. A radial profile of the excess counts per square degree 
in rings of 0.1° width starting from the Westerlund 1 position is shown in Fig. 
SI The VHE 7-ray emission from within the optical boundary of the stellar 
cluster (the innermost ring) is consistent with the overall picture of a rather flat 
emission out to a distance of ~ 0.9° from the Wd 1 position (a fit of a constant 
yields a moderate x 2 / n df of 14.6/8). Fig. [3l right shows preliminary H.E.S.S. 
significance contours from 4 to 8 a in steps of 1 a after integrating events within 
0.22° on top of the HI channel ma p at a velocity of -55 km s~\ corresponding to 
the distance to Wd 1 as derived bv lKothes Doughert y (2007). The comparison 
may suggest a correlation between the VHE 7-ray emission and the HI data in 
some parts. Further spectral and morphological studies of the H.E.S.S. emission 
and observational data in radio and X-rays could help to shed some light on the 
origin of the TeV 7-ray emission and its possible connection to the super star 
cluster Westerlund 1. 



6. Discussion 

Massive stellar clusters provide sufficient kinetic energy in stellar winds, col- 
liding wind binary systems and from collective wind and/or wind/SNe ejecta 
effects to accelerate particles to relativistic energies. After the detection of VHE 
7-rays from the direction of the open stellar cluster Cyg-OB2 by the HEGRA 
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Figure 4. Preliminary radial profile of the uncorrelated 7-ray excess sky 
map starting at the Westerhmd 1 position (white star in Fig. [3j left). 



Collaboration ( Aharonian et al. 2005), H.E.S.S. discovered a source, coincident 



with the stellar cluster Westerund 2. Furthermore, one of the brightest uniden- 
tified H.E.S.S. sources, HESS J1614-518, seems also be connected to a stellar 
cluster, Pismis 22. Westerlund 1 as the most massive stellar cluster known in our 
galaxy is a perfect target for VHE observations given its massive star content, 
age and binary fraction. The detection of extended emission from the vicinitiy 
of Westerlund 1 has been reported in this work. All these results suggest that 
multi-TeV particle acceleration may be linked to several massive stellar clusters. 
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